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Abstract 

Fusion can be described by the time evolution of a dinuclear system with two degrees of freedom, 
the relative motion and transfer of nucleons. In the presence of the coupling between two collective 
modes, we solve the Fokker-Planck equation in a locally harmonic approximation. The potential 
of a dinuclear system has the quasifission barrier and the inner fusion barrier, and the escape rates 
can be calculated by the Kramers’ model. To estimate the fusion probability, we calculate the 
quasifission rate and the fusion rate. We investigate the coupling effects on the fusion probability 
and the cross section of evaporation residue. 
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I. INTRODUCTION 


Fusion reaction has been greatly investigated, and superheavy elements have been syn¬ 
thesized in experiments. Many theoretical models describe heavy ion reaction and estimate 
the cross section of evaporation residue. Although much effort has been devoted to un¬ 
derstanding fusion, it is still challenging to analyze and to predict data quantitatively!!]. 
Especially, the fusion probability is the least understood factor contributing to the cross sec¬ 
tion. The fusion probability has been studied as diffusion, and it is expected to be a logistic 
function[2H3]. However, it might be sensitive to details of dynamics and be important for 
more quantitative analysis of data. In this work, we investigate the fusion probability by 
using the Fokker-Planck equation in the dinuclear system concept. 

A dinuclear system consists of two colliding nuclei captured in an almost touching 
con£guration[5|. In this concept, fusion can be understood as the time evolution of a 
dinuclear system. This description is different from other models in that it retains the 
individuality of nuclei and considers quasihssion. A dinuclear system has two degrees of 
freedom, the relative motion and transfer of nucleons. To describe the relative motion, the 
internuclear distance R is used as a collective coordinate. For transfer of nucleons, we define 
the mass asymmetry parameter 


Ai — A2 
Ai A2 


( 1 . 1 ) 


where Ai and A 2 are mass numbers of two nuclei. Each mode of motion is a diffusion process 
in a potential with an energy barrier. The escape rate over the barrier can be calculated by 
the Kramers’ model. There has been much research on the Kramers’ model (for a review 
and references, see Ref. 0 )- The extension of the model can be applied to fusion reaction. 

The time evolution of a dinuclear system is described by the Fokker-Planck equation. For 
asymmetric reactions with not too high mass asymmetry, we have the following equation |21 

m: 


dl 
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( 1 . 2 ) 


where /ij, 7 *, and Di {i = R, p) are the diagonal components of the inverse mass tensor, the 
friction tensor, and the diffusion tensor, respectively. By the fluctuation-dissipation theorem, 
Di = PiTi, where T* = {hui/2) coth.{hui/2T) is the effective temperature with the energy of 
zero oscillations fkiJi/2 and the thermodynamic temperature T = ^JE*/a (a = A/12 Me 
is the level density parameter). 

The Kramers’ problem described by Eq. (1.2) is useful in other quantum mechanical 
systems as well as a dinuclear system. The escape rate in two-dimensional potential has 
been calculated^ E], but we use a different strategy to calculate the diffusion rates of two 
collective modes. By assuming that two modes are weakly coupled through the inverse 
mass tensors, we solve the Fokker-Planck equation perturbatively. We then calculate the 
quasihssion rate and the fusion rate to estimate the fusion probability. Although we apply 
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FIG. 1. A schematic potential for the Kramers’ model. For moderate to high friction, a system 
at a local minimum x = Xa can escape the energy barrier (Eb > T) at x = Xb to go to the lower 
bound state at x = Xc- 


the rates to quasifission and fusion, our results for the diffusion rates can be applied to any 
escape problem in the presence of friction and diffusion. 


In Section we review the Kramers’ model to calculate the escape rate. In Section 


we extend the model with two collective modes which are coupled through the inverse mass 
tensors. For weak coupling, we solve Eq. (1.2) perturbatively and calculate the diffusion 
rates of two modes. By applying the escape problem to a dinuclear system, we consider 
asymmetric fusion reactions in Section IV We review the dinuclear system concept and 
estimate the coupling effects on the fusion probability and the cross section of evaporation 
residue. Finally, we summarize our results in Section |Vj 


II. KRAMERS’ MODEL 


The fusion probability can be estimated by the Kramers’ model. Kramers considered 
the escape problem as a one-dimensional Brownian motion in a potential with deformation 
energy [To]. In this section, we review the Kramers’ model, and the extension with two 
degrees of freedom will be discussed in the next section. 

We consider a system initially at a local minimum x = Xa (see Fig. [^. For moderate to 
high friction, the system can be excited enough to overcome the energy barrier at x = Xb 
and go to the lower bound state at a; = Xc- Since the energy barrier is supposed to be higher 
than the thermal energy {Eb > T), the escape process is slow. It can be understood as 
quasistationary diffusion. 

The evolution of the system is described by the Fokker-Planck equation[TU] 


d d dU{x) d 

dt dx ^ dx dv 


f{t,x,v) 


d 


vf{t,x,v) + fiT 


df{t,x,v) 

dv 


( 2 . 1 ) 


In the diffusion limit, the system is in equilibrium everywhere except near the barrier at x = 
Xb- Thus, the dynamics may be evaluated by considering the potential around this point. In 
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a locally harmonic approximation, the potential is given by U{x) = —ojI{x — XbY/^n + U{xh), 
where ujb is the angular frequency of the unstable state at the barrier. Around the barrier, 
there are neither sources nor sinks, and the distribution function satishes the stationary 
Fokker-Planck equation. 

In the stationary limit, we solve the Fokker-Planck equation by using a Maxwellian dis¬ 
tribution 


/(x,u) = C{.x,v) exp 



+ U{x) 



( 2 . 2 ) 


With the boundary conditions v) 

isP [TU] 


C(a:,u) 




1 and C{xc, v) = 0, the Kramers’ stationary solution 

(2.3) 


POO 


/ du exp 

' {x-Xb)-\v/ul 

2^T'yX 


where A = ( 7 / 2 )^ + — 7/2. The stationary diffusion current is calculated by 


J(xb)= / dvvf{xb,v), 

J —CO 

and the number of particles at the initial state is 




N'{xb)= / dx dv f{x,v), 


' —CO J —CO 


(2.4) 


(2.5) 


where we use the locally harmonic potential U{x) = ul(x — a:a)^/2/i -|- U(xa) around x = Xa- 
The escape rate is then given by the flux over population [TO] 


■pfcr _ i^b) 

~ M{xb) 


XcJn 


2nuJb 


exp[-Eb/T] , 


( 2 . 6 ) 


where the energy barrier is Eb = U{xb) — U{xa)- 

The Kramers’ rate with the parabolic potential might give an upper bound of the es¬ 
cape rate mm- For high friction, there can be recrossing over the energy barrier, and 
anharmonic corrections of the potential barrier reduce the escape rate. In Ref. [ 12 ], the lead¬ 
ing order effects of anharmonic potentials have been calculated perturbatively. In the next 
section, we use the same strategy to solve Eq. (1.2) in a locally harmonic approximation. 


III. COUPLING THROUGH INVERSE MASS TENSORS 


The Kramers’ quasistationary rate can be applied to two-dimensional potential of a dinu- 
clear system. In general, the radial mode and the mass asymmetry mode are not separable, 
and they are coupled with each other. The coupling between two collective modes is weak 
for almost symmetric reactions, but it grows as the asymmetry increases [T3|. For simplicity, 
we consider asymmetric dinuclear systems where the coupling between two modes is still 
weak. By following Ref. [T2] in which the Fokker-Planck equation is solved perturbatively, 
we investigate the coupling effects on the escape rates. 


Similar to Eq. (2.2), we consider a distribution function of the form 


f{R,r],VR,Vrj) = C{R,r],VR,Vrj) exp [-{v\/2hrTr + vy2iJrjTrj + U{R,r])/Teff)] . (3.1) 
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Around the barriers {R = Rb, V = Vb), the two-dimensional potential is assumed to be 


U{R-,v)/Teff — —u'ji^{R — Rb)‘^/2fijiTji + Ueff{Rb)/Tji 

~ ~ Vb) + Ueff{rib)/Tri , (3.2) 

where Ueff{Rb) and Ueff{rib) are the effectively one-dimensional potentials in each collective 
coordinate. If there is no coupling through the inverse mass tensors, the escape rate is given 
by the Kramers’ rate in Eq. (2.6) for each degree of freedom. In that case, we have the 

(3.3) 


zeroth order equation of motion 


CoCo{R,V,vr,v^) = 0 , 


where we have defined 


^ d dU d ^ d *rT. 


dvl 

d dU d ^ d ^ 


with 7 * = pj 7 j {i = i?, 7 ). Since R and 7 modes are separable at the barriers, we have 

Co(fl, i>,) = C,R(R,VR)C,„(ri,v^), 


(3.5) 


and each factor is given by the Kramers’ stationary solution in Eq. (2.3): 

,2 poo 


Cr{R:Vji) — 


Cviv,Vn) = 


Ul 


Rb 


\/27r/r/jTR7^AK J{R-Rf,)-\jiVji/ujj^ 


du exp 


UJ% 

Kb 


uj: 


Vb 


y^2Tl^j^Tri'y*Xri J (^ri-rii,)-\^Vr,/uj?i^ 


du exp 


‘2fJ.RTR'y^\R 




_ ‘2lJ'r}Tri'y^\ri 


(3.6) 


where A, = ^^^72)^+^ - 7,72 (* = R,v). 


When the relative motion and the mass asymmetry are coupled by Eq. (1.2), we define 
the first order operator 


A = 


dfiR 


drj Hr 


1 1 dUr] f^R 2 

VRVr^ + V. 


2 dR ^ 




d 


dv 


R 


dflr] 1 


VRVr^ + 


1 dpR Hr, 


dRni"'''^ ' 2 dT] hr^ 

By expanding the solution in the order of coupling. 


A 

dVr, 


(dHr, 1 
\dR Hr, 


Vr + 


dHR 1 \ 

dv X^R^V 


(3.7) 


C{R,V,VR,Vr,) = Co{R,V,VR,Vr,) + Cl{R, V, VR, Vr,) + ( 2 (7?, ??WRWr,) H-, (3.8) 

we have the equations of motion [ 12 ]: 

- /:oCi(7?,7WRWr,) = RiCoiR^V^VR^Vr ,), (3.9) 

-£ 0 ( 2 ( 7 ?, 7 , Vr, Vr,) = £lCl(7?, 7, Vr, Vr,) , (3.10) 


5 



























and similarly for higher orders. In the following snbsections, we will solve the first order 
eqnation of motion. 

After solving for (i{R,ri,vji,Vrj), we need a strategy to calcnlate the escape rate in the 
two-dimensional potential. By noting that the zeroth order resnlt shonld be two independent 
escape rates, we calcnlate the diffnsion cnrrent and the nnmber of particles as follows. For 
the relative motion, the current at the barrier is 


JhiRb-, Vb) — 


f*oo poo 

dv] 


dVr, VRf{Rb,r]b,VR,v ^), 


jR I iMurf 'jbi '^ni '^r/J 

' —OO J —CO 


and the number of particles is dehned as 


Rl'ni.Rb.Vb) — 


"Rb 


dR 


dvf 


JR I J y^'^i 'Ibi '-'Hi '-"q) 

'—OO J—OO J—OO 


‘il'n f{R< 1b, Vk, vJ , 


(3.11) 


(3.12) 


where we have set r] — rj^to cancel eiq)[—U„ff{lb)/T,i] when taking the flux over populatioiU 
Similarly for the mass asymmetry mode, we have 


Vb) 


f *00 poo 

dv] 


dVri Vrjf{Rb,r]b,VR,V,^) , 


jR I U'qj Uyj) 

' —OO J —OO 


and 


/ OO nrny noo 

dvR dr] dVrjf{Rb,T],RR,Rv)- 

■OO J —OO J —OO 


Then the escape rates of two collective modes are given by the flux over population; 

7 (-^^7 


p jRiRbi Vb) 1 

^R = TrTR -f 

J^RyRbi Vb) 


Afr^{Rb,Vb) 


(3.13) 

(3.14) 

(3.15) 


Depending on reactions, not every term in Eq. (3.7) might be important to calculate the 


diffusion rates. To investigate the effect of each coupling, we solve the hrst order equation 


of motion separately. We present details in Section III A and only the results are presented 
in the following subsections. 


A. Cl — 


IS 


In this subsection, we assume Ci = Since the right hand side of Eq. (3.9) 

(3.16) 


RlCoiRlVl'URl'Uq) = -- CviVRq) 1 


we 


expect the hrst order solution to bfi 


3l 


R 


Ci{R,v,VR,Vq) = 


d]iR 1 
dv V'R 


p/D ldCR{R,VR) 

Pl{Ri Vi VRi Vrj) -- CqiVi Vq) 


R 


+ P2{RiVRRiVq) 


dCRiPi Vr) dCqiViVq) 


dv 


R 


dv„ 


, (3.17) 


^ Due to the special way to calculate the number of particles with two degrees of freedom, we calculate the 
number by using only the zeroth order solution. 

^ In Eq. (3.17), there is a trivial term a(R{R,VR) (^r]{Vii’q)i where a constant is chosen to be a = 0. 
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where Pi{R,r],VR,Vr^) {i = 1, 2) is a polynomial of collective coordinates and velocities. By 
plngging Eq. (3.17) into Eq. (3.9), we have 


'’’‘Wr ■ 


8 = 




-'R 


d d d 


dVn 


j 


and 


d d d 

+ '1*r)^R 


Pi{R, r], vr, v^) = -VRVr, , (3.18) 




R 


dvR 


d d d d‘^ 


'dvr, 


dv^ 




P2{R,r],VR,v^ 
dPi{R,r],VR,Vr,'^ 


dvr. 


. (3.19) 


The polynomials are fonnd to be 

Pi{R,ri,VR,Vn) =ai{R- Rb){rj - r]b) + a 2 {R- Rb)vr^ + 03(77 - r]b)vR + a^VRVr ^, 
P2{R,ri,VR,Vn) = a^{R - Rb) + RqVr , (3.20) 

where a* (i = 1, • • • , 6 ) is a constant. It is straightforward to determine the constants: 

^lb^‘Rb^'y*r) + 7r + 4AR)a4 


Oi = 


02 — 


03 — 


04 = 


Os — 


06 — 


A_r( 3A^ — 4A|.) + 7i)[7r?(-^?? ~ 2Ai?) + XR^SXr/ — 4Ai{)] + 7^(A^ — 2A^) 
_ ^Rt^R^R + 3Afi + + ujjiJaA _ 

7h( 2A|. — A^) + A/j(4A^ — 3A^) + 7i)[7fi(2AR — A,;) + Aij(4A/j — 3Ar,)] 

_+ %i\ ~ -^r) + 7r-^r]Q4_ 

A^(3A2 - AXl) + 7;[7le(A, - 2Xr) + Xr{3X, - 4A^)] + ^r^ - 2X1) ’ 
A^Aj;[Afi(3A^ - 4A|) + 7)){7R(Ar; - 2Xr) + Ar(3A,, - 4Ah)} + 7fl(A^ - 2A|)] 
(A^ - 2Xr){2Xr - A^ + 7|j)(2A^Ar + uj'^J[Xruj'^^ + A,,(A| + co^J] 
27j)/^»7^r;[(7R + A,; + 3 Ah)02 + (^R^Ra] 


A^ + 4:X^Xr + 3A^ + 7^(A^ + A/j) + ’ 

2'y^fJ-riTrj[—a2 + (A^ + A/j)o4] 

(A,, + Ai?)( 7 |j + A^ + 3 Ak) + uJr^ 


(3.21) 


By using Eqs.(|3.11|) and (|3.12|), we find the diffusion current and the number of particles 

d/RR 1 


for the relative motion to be 


jR{.Rh^ Vb) — 


TT^riTri ^rTrXr 


^Rb 


1 + 


dv l^R 


(o4/ijjX)y T Ojq) 


_2_ 1*r^r\ 

nUnTri 


X exp [-Ueff{Rb)/TR - Ueff{vb)/Tr,] , 

N'R{Rb,r]b) = \/2TTIJ,nTrj ^^^'^^ exp [-Ueff{Ra)/TR - Ueff{r]b)/Trj] . 


^Ra 


(3.22) 
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Then the flux over population is used to determine the escape rate 

1 




1 + 


dv fJ'R 


{.CLil-^rjTrj T ^ 0 ) 




(3.23) 


where we have used the Kramers’ escape rate in Eq. (2.6). Similarly for the mass asymmetry 
mode, we have 


Vb) 




U). 


exp [-U^ff{Rh)/TR - UeffiVb)/Trj] , 


Vb 


Af,{Rb, Vb) = V27rpRT^^^ exp [-Ueff{Rb)/TR - U,ffiva)/T,] , (3.24) 

and the diffusion rate is 




Va 


F — F* 

J- r? ^ n 


(3.25) 

We note that the escape rate of the R mode changes by the coupling, while the rj mode is not 
affected. The coupling effects depend on the friction coefficients and the angular frequencies 
at the barriers. 


p> r _ 2 _^ 

- 2 dRR^^Vdvn 


The first order solution is of the form 


Ci(i?,?7,ni?,n^) = - 


1 O^'ri V'R 


2 dR 


p/D ^dCRiR,VR) 

Pi {R, p, Vr, Vr,) - — - Cviv, 


R 


+ P2{R,V,VR,Vrj) 


OCr^R, vr) dCrj{v,'Vn 


dv 


R 


dv„ 


and the polynomials are found to be 

Pi{R,r],VR,Vr^) =bi + b2{r]- Vb)^ + bsiv - Vb)vr, + 
P2{R,r],VR,Vr^) =b^{v- Vb) + baVr ^, 
where the constants are 

2Hrj'y*Tr,bi 


bi = 

b2 = 
bs = 
64 = 


A 


R 


2ut bi 


Vb 


(A 2 j- 2 A 2 )+ 7 ;(Ah- 2 A,)’ 

“^XROO^^bi 

-(Xl-2Xl)+r,{XR-2X,)^ 
(A| - 2A^) + 7*(A/j - 2Xn) 


{% + Xr){Xr — 2Xri){2'y* + A/j + 2Xrj) 

^ _ 2p^7*T^[(7* + Ar + 3A,,)63 + 2a;^^&4] 

X% + 4ArA^ + 3A^ + 7 j)(AR + Xrj) + u, 
2/^»?7j)^r;[~&3 + 2 (Ar + A,,) 64 ] 


2 ’ 
Vb 


bf, = 


{Xr + A,,) ( 7 * + Ar + 3A,,) + 


(3.26) 


(3.27) 


(3.28) 























The escape rates of two collective modes are 


Tr = r 


kr 


r = r 

-L 7? -L r; 


kr 


1 + h + W, ( 2 + 




0J„ 


+ be 


uJi 


dR 2 tiTrUJr^ ^ X / -rj. 

In this case, the diffusion rate of the rj mode is affected by the coupling. 


(3.29) 

(3.30) 


n _ djin 1 d 


This case is similar to Section III A by interchanging {R,vr) and (p,n^). The first order 
solution has the form 


Ci{R,ri,VR,v^) = 


djXn 1 
dR firi 


Pi{R,r],VR,Vrj)CR{R,VR) 


dCr^{ri,Vri 


dVr. 


+ P2{R,r],VR,Vr, 


SCr^R, Vr) dCr,{r],Vr, 


dv 


R 


dVr, 


, (3.31) 


where the polynomials are 


Pi{R,r],VR,Vr,) = ci{R - Rb){r] - r]b) + C2{r] - r]b)vR + C3{R - Rb)vr^ + C4VRVr^, 

P2{R, r], Vr, Vr,) = Cb{r]- r]b) + c^Vr, . (3.32) 


The constant q (i = 1, • ■ ■ , 6) is given by Eq. (3.21) except that (jr, Xr, ojr^, ^r,Tr) and 
(7*, Xr,, uJr,^, iJir,,Tr,) interchange. The escape rates of two collective modes are 




kr 
V ~ ri 


r„ = r 


1 + 


d/Xr, 1 

dR fir. 


{cifiRTR + ce) 


2 ArA^ 7* 

T^fXRTR UJR^U}^^ 


(3.33) 

(3.34) 


D. 


C 


_ 1 d^iR biri 2 d 

1 2 dr, ^R dvr, 


This case is similar to Section IIIB by interchanging {R,vr) and {ri,Vr,) 
solution is of the form 


The hrst order 


Ci{R,ri,VR,Vr,) = - 


f dfiR fir, 

2 drf fi\ _ 


Pl{R,rf,VR,Vr,)CR{R,VR) 


dCn{rf,Vr, 


dvn 


+ P2{R,rf,VR,Vr,) 


dCR^P, Vr) dCr,{rf,Vr, 


dv 


R 


dVr, 


(3.35) 


with the polynomials 


hi Vr, Vr,) — di + d2{R — PbY + '^3(-R ~ Rb)vR + d^V^ , 

P 2 {R, rf, Vr, Vr,) = de{R - Rb) + deVR . (3.36) 
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The constant (ij (i = 1, • • ■ , 6) is given by Eq. (3.28) except that (7 ^, Ar, Ci;rj,pr,Tr) and 
(7*, A^, interchange. The escape rates of two collective modes are 


rR = r 


kr 


1 + 


d^R 1 
dri /i^' 




+ «6 


CJ 




u 


Rb 


T = T 

J- r; J- r; 


kr 


, (3.37) 
(3.38) 


F JLii„ nr- — QAtfi 1 

hj. 1-1 — Ur or 1-1 — 




We assume Ci = ~'^^'Vr- The first order solution is 

d 1 

Ci(-R,7WrW7?) = -^—PiR,V,'i^R,V'n)(Ri^^VR)C^{ri,Vr^), (3.39) 

Cy-t\ 

where the polynomial is found to be P{R, rj, vr, u ^) = —(i? — Rb) + c. For a constant c = 0, 
this coupling does not change the diffusion rates of two modes. 

Similarly, the coupling Ci = does not affect the escape rates. 


F. Coupling Effects on Escape Rates 


As seen in Eqs. (3.23), (3.30), (3.34), and (3.37), the escape rates are affected by the 


coupling between two collective modes. The coupling effects generally depend on the friction 
coefficients (7R,7j)) and the angular frequencies (u;Rj,a;^J. In this subsection, we investigate 
the degree of the coupling effects on the escape rates, by using specihc parameters. For a 
dinuclear system, the diffusion rates of the R and rj modes are interpreted as the quasihssion 


rate and the fusion rate, respectively (see Section IV). 

The inverse mass tensors of a dinuclear system are given 


^^R 


Am 1 — 7]"^ 

1 z/ 

Am 2^/^b^^ ’ 


1 - 


u 


1 — 7 ^ 


with ^(^o-^ih^)(l-^s), 


(3.40) 


where m is the nucleon mass, s = R — Ri — R 2 , ^0 = 16, .^1 = 17.5, ^ = 0.3 fm“ , and 
b = 1 fm. For a dinuclear system with A = 300, we use the coefficients of hwR^ = 2MeV, 
huj^ji = 1 MeV, and = hy* = 2 MeV. 

In Fig. we present the ratios of the escape rate to the Kramers’ rate. The blue solid 


the fusion rate of Eq. (3.30) in Section IIIB 


lines are for the quasihssion rate of Eq. (3.23) in Section III A, the violet dashed lines are for 


the red dashed-dotted lines are for the fusion 


rate of Eq. (3.34) in Section III C| and the green dotted lines are for the quasihssion rate of 
Eq. (3.37) in Section HID Fig.(a) shows the coupling ehects depending on the excitation 
energy for s = 1.5 fm and rj = 0.6. For the particular choice of parameters, the fusion rate 


in Section IIIB decreases by 10-25%. On the other hand, the coupling in Section IIIC 


increases the fusion rate by 20-25%. Thus, in the presence of both couplings, there might 
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FIG. 2. The ratios of the escape rate to the Kramers’ rate (a) as a function of the excitation energy 
for s = 1.5 fm and rj = 0.6, (b) as a function of the distance for E* = 30 MeV and r] = 0.6, and 
(c) as a function of the mass asymmetry parameter for E* = 30 MeV and s = 1.5 fm. The blue 
solid lines are for the quasifission rate of Eq. (3.23) in Section III A, the violet dashed lines are for 
the fusion rate of Eq. (3.30) in Section IIIB, the red dashed-dotted lines are for the fusion rate of 
Eq. (3.34) in Section IIIC and the green dotted lines are for the quasihssion rate of Eq. (3.37) in 
The parameters of = 2 MeV, = 1 MeV, = h'y* = 2 MeV, and A = 300 


Section 


HID 


are used. 


be cancellation in some of the effects. For the quasihssion rate, the coupling effects are 
relatively low. The coupling in Section |IIID increases the rate by 2.5%, and the effect in 
Section III A is less than 1%. Fig. (b) shows the coupling effects depending on the distance 
for E* = 30 MeV and r] = 0.6. As the distance increases, the effects on the fusion rate grow 
due to the increasing mass of the rj mode. The fusion rates change by 10-30% while the 
quasihssion rates are not ahected much. The couplings ehects also depend on the mass 
asymmetry as in Fig. (c). Dinuclear systems with higher mass asymmetry are ahected 
more by the coupling through the inverse mass tensors. 
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IV. DINUCLEAR SYSTEM 


In this section, we review the dinuclear system concept by following Refs. laEHn] and 
apply the Kramers’ escape problem to a dinuclear system. We discuss the diabatic nucleus- 
nucleus potential and the driving potential. By the Kramers’ model, quasihssion and fusion 
are described as diffusion in two collective coordinates. The cross section of evaporation 
residue is determined by estimating the fusion probability and the survival probability. 
Based on the study in the previous section, we discuss the coupling effects on the fusion 
probability and the cross section. 

A. Diabatic Potential 

The diabatic nucleus-nucleus potential consists of the Coulomb potential, the nuclear 
potential, and the centrifugal potential: 


V{R, J) = Vc{R) + Vn{R) + VrotiR, J ), 


(4.1) 


where J is the angular momentum. We consider axially symmetric systems described by 
Ri{6i) = i?oi [1 + A^ 2 o( 6 'i)], where (3i is the quadrupole deformation parameter and Roi = 


I/O 

roA/ with ro = 1.16 fm. Then the Coulomb potential is given bv[T 8 ] 



(4.2) 


The centrifugal potential is 



(4.3) 


where X is the moment of inertia. 


The nuclear potential in the double folding form is 



(4.4) 


where F{ri — r 2 ) is the nucleon-nucleon interaction. A well-known ansatz for the density 
dependent interaction gives [T9] 



(4.5) 


where 



(4.6) 


and similarly for F^x- The density p{r) = pi{r) + Pii^) is given by the Wood-Saxon form 



(4.7) 
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R T| 

FIG. 3. (a) The diabatic nucleus-nucleus potential has an energy barrier, = V{Rb) — V{Rm)- 
(b) The driving potential has the inner fusion barrier, Bj^^ = U{r]BG) ~ U{rii). There is another 
barrier for symmetrization of a dinuclear system, Bgym = [/(r/i+i) — U{r]i). The quasifission barrier 
is defined as Bgj = Min[i?^, Bsym]- 

We use the parameters of Cq = 300 MeV fm^, = 0.09, = 0.42, = —2.59, = 0.54, 

Po = 0.17fm“^, and oq = 0.55fm. 

For J = 0, a schematic nucleus-nucleus potential is shown in Fig. (a). In contrast to 
the adiabatic potential, the diabatic nucleus-nucleus potential diverges at short distances. 
It has an energy barrier, = V{Rb) — V{Rm)- The potential energy of a dinuclear system 
is 

U{R, f], J)=Bi + B2 + V{R, J) , (4.8) 

where Bi {i = 1, 2) is the binding energy including shell effects. A schematic of the driving 
potential is shown in Fig. (b). It has the inner fusion barrier, Bj^^ = U{r]BG ) - U{m), 
which is the difference between the driving potential at the Businaro-Gallone point rj = rjBG 
and the initial value at r/ = rji. B'^y^ = — U{r]i) is the barrier for symmetrization of 

a dinuclear system. The quasihssion barrier is dehned as Bgf = Min[i?^, B]y^] 

In the dinuclear system concept, fusion reaction can be described as follows. The initial 
dinuclear system is supposed to be at the bottom of the potential pocket since the system at 
the minimum survives longest. If the system overcomes the quasihssion barrier, it decays. 
The escape over the barrier can be understood as diffusion in the relative distance. On the 
other hand, if a dinuclear system overcomes the inner fusion barrier, it undergoes fusion and 
forms the compound nucleus. Fusion can be understood as diffusion in the mass asymmetry 
parameter through transfer of nucleons. There is competition between quasihssion and 
fusion. The fusion probability can be estimated by the Kramers’ model. After fusion, 
the system is highly excited and is likely to hssion. However, there is a possibility that a 
dinuclear system remains as evaporation residue by emitting neutrons. This possibility is 
dehned as the survival probability. In the following subsections, we discuss how to calculate 


Unlike the assumption made in Section 


III 


R and rj modes are coupled with each other through the 


potential even if there is no coupling through the inverse mass tensors. 
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the cross section of evaporation residue by estimating the fusion probability and the survival 
probability. 


B. Evaporation Residue 


The cross section of evaporation residue is determined by the effective capture cross 
section, the fusion probability, and the survival probability [20]: 

P(j^(yE 

cm )Wsur{E,m): (4.9) 

with 


O', 


cap 




vrh" 


2 /rEe 


-{Jmax + E)‘^E{Ecm) , 


(4.10) 


where /r is the reduced mass, Ecm is the bombarding energy in the center of mass system, 
and E{Ecm) is the transmission probability. 

The fusion probability describes the possibility for a dinuclear system to overcome the 
inner fusion barrier, competing with quasihssion. We apply the Kramers’ escape problem to 
a dinuclear system with the quasihssion barrier, Bqf = Ueff{Rb) — Ueff{Ra), and the inner 
fusion barrier, Bf^s = UeffiVb) — Ueff{j]a)- Then the fusion probability can be estimated 

aflini [m [21] 

= ( 4 - 11 ) 


+ Tr; 

the coupling effects on the escape rates depend on the friction 
coefficients and the angular frequencies at the barriers. For the particular choice of param- 


As discussed in Section III 


eters in Section III F , the coupling in Section III B reduces the fusion rate, so we expect 
the fusion probability to decrease. On the other hand, the coupling in Section III C might 
increase the fusion probability. The effects on the quasihssion rate in Section III A and 
Section ^ might be ignored, in comparison to the fusion rate. 


Although the values of the coefficients are not known, a phenomenological formula for 
the fusion probability has been obtained [20]: 

p 1.25 exp - Bqf)/TDNs\ . . 

l + 1.25exp[-(5;„,-i?,,)/T^^,]’ ^ 

where the temperature of a dinuclear system is Tdjvs = a/ E*^j^g/a with the excitation 


energy, ~ UiVi)- Based on the study in Section III, the coupling ehects might 

be included in the coefficient of the exponential term in Eq. (4.12). Thus, the variation of 


this coefficient allows us to estimate the degree of the coupling ehects. 

After fusion, the compound nucleus can remain as evaporation residue by emitting neu¬ 
trons, not to hssion. This possibility is described by the survival probability. The survival 
probability can be estimated to b(0221 123] 

^n{Ecj^i) 


IF. 


'(Ecn) — Pxn{E'^j^) 


=1 ^/(-^CV,i) + 


(4.13) 


Fij -|- is the sum of the diffusion rates of R and 77 including symmetrization of a dinuclear system. 
Except fission and neutron emission, we ignore all other processes such as a particle emission and 7 ray 
emission. 
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where P^niE^j^) is the probability of emitting x neutrons at the excitation energy E^j^, 
ris the width of hssion at the excitation energy E^j^ and Vn{E^^^) is the width 
of neutron emission at Since an emitted neutron carries the average energy of -B„ + 2T 

{Bn is the neutron separation energy), we have ~ {Bn,i + 2Tj), where 

77'* _ TTi* 

^CN,l — ^CN' 

For X >2, the probability of x neutron evaporation is 


P^n{EcN) = P[x] - P[X + 1] , 


where 

P\x] = 1 


2x-3 


,-A./r 


i + E 




with 


2=1 




— Eqj^ 


(4.14) 


En,i ) (4-15) 

2=1 


and the effective temperature is T = E^j^/1.5 a. For In channel, we use 

-Pln(-E'C'w) ~ [~{EcN ~ En ~ 2T) /2(T ] , 


(4.16) 


where T = and a = 2.5 MeV. The ratio of the neutron emission width to the 

hssion width is given hymm 


\ ^ AA^/\E*c^ - Bn) 

f fc [2 — 1] 


exp 


2-\/a ( ^ V^CN ^ 


(4.17) 


Here, k = 9.8 MeV, and the hssion barrier depends on the excitation energy as Bf{E^j^) = 
Bf{EcN = 0) exp [-E^N/Ed] with Ed = 0.4 H^/3/a|25]. 


C. Coupling Effects 


In this section, we investigate the degree of the coupling ehects on the fusion probability 
and the cross section for the reaction ^gFe+^glPb —Hs + xn {x = 1, 2). We use 9i = 0, 
02 = vr, Jmax = 10, and T{Ecm) = 1 for estimation. The nuclear data such as the Q value, 
the hssion barrier, and the neutron separation energy is from Refs. [2^ 127] . 

Fig-i a) shows the fusion probability depending on the excitation energy. The solid line 
is calculated with the phenomenological formula in Eq. (4.12), the dashed lines include 10% 


coupling ehects compared to the phenomenological formula, and the dotted lines include 
20% coupling ehects. We mention that 10-20% ehects are caused by corrections to both 
the quasihssion rate and the fusion rate. The fusion probability generally grows as the 
excitation energy increases. The degree of increase depends on the degree of couplings. 
Fig. I^b) shows the cross section of evaporation residue calculated by Eq. (4.9). The red 
lines are for the In evaporation channel, and the green lines are for the 2n channel. The solid 
lines are calculated with the fusion probability given by the phenomenological formula. The 
dashed and dotted lines are calculated with the fusion probability including 10% and 20% 
coupling ehects, respectively. As the excitation energy increases, the cross section initially 
grows and then decreases due to the survival probability. At the maximum of the cross 
section, the coupling ehect is largest. The black points show the experimental data with 
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FIG. 4. (a) The fusion probability for the reaction ^gFe Pb 


Hs -|- xn. The solid line 


represents the phenomenological formula given by Eq. (4.12), the dashed lines include 10% coupling 


effects compared to the phenomenological formula, and the dotted lines include 20% coupling 
effects, (b) The cross section of the evaporation residue for the reaction IgFe-l-^glPb —Hs-|-xn. 
The red lines are for the In channel, and the green lines are for the 2n channel. The solid lines are 


calculated with the fusion probability given by Eq. (4.12), the dashed lines include 10% coupling 


effects on the fusion probability, and the dotted lines include 20% coupling effects. The black 
points represent the experimental data with errors[2T]. 


errors izj. The cross section calculated with the phenomenological fusion probability agrees 
with the experimental data within the errors. However, the coupling effects might produce 
uncertainties on the cross section. We note that the magnitude of the error is comparable 
to 10 - 20 % coupling effects at the maximum of the cross section. 

We should mention that there are other factors contributing to the uncertainties of the 
cross section, in addition to the coupling effects between two collective modes of a dinuclear 
system. Especially, the degree of the survival probability highly depends on calculation 
methods and parameter values. However, in Fig. we present the degree of the coupling 
effects for a reaction in which the calculated cross section agrees with the experimental data. 


V. SUMMARY AND DISCUSSIONS 


In this work, we considered the Kramers’ escape problem with two collective modes and 
calculated the diffusion rates in the weak coupling limit. Due to the coupling through the 
inverse mass tensors, the rates have corrections depending on the friction coefficients and 


the angular frequencies at the barriers. Our results are given by Eqs.(3.23), (3.30), (3.34), 


and (3.37). They can be generally applied to any escape problem in a quantum mechanical 


system with friction and diffusion. 

We have applied the escape problem to a dinuclear system to estimate the fusion prob¬ 
ability. For the coefficients of = 2 MeV, hwn^ = 1 MeV, and = hj* = 2 MeV, we 
present the coupling effects on the diffusion rates in Fig. The fusion rate changes up to 
30% while the quasihssion rate is barely affected by comparison. Since the fusion probability 
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is determined by the diffusion rates of two collective modes, the couplings affect the fusion 
probability and the cross section of evaporation residue. Fig. shows 10-20% coupling 
effects for the reaction ^gFe Pb —Hs + xn {x = 1, 2). The effects depend on the 
excitation energy, and they are the largest at the maximum of the cross section. While the 
calculated cross section agrees with data, the experimental uncertainties are comparable to 
10 - 20 % coupling effects at the maximum of the cross section. 

We should mention that there are other uncertainties in estimating the cross section of 
evaporation residue. Some of the effects might be greater than the coupling effects on the 
fusion probability. Especially, the survival probability highly depends on the calculation 
methods and parameter values. However, the purpose of this work is to investigate the cou¬ 
pling effects on the fusion probability for more quantitative analysis of data. We calculated 
the leading order corrections to the quasifission rate and the fusion rate. The effects of each 
coupling might be important to understand fusion depending on reactions. Since we do 
not have complete understanding of fusion dynamics nor nuclear data in details, we need a 
phenomenological strategy to compare the analytical results with experiments. 

We focused on asymmetric dinuclear systems in the presence of the coupling through the 
inverse mass tensors. We need many assumptions to obtain Eq. (1.2), but the general de¬ 
scription of a dinuclear system involves other effects. First, we ignored dfiji/dR and dfi^/drj 
which might be important in some reactions. Second, we have non-diagonal components of 
the inverse mass tensor, the friction tensor, and the diffusion tensor as well as the diagonal 
components of the kinetic energy. By including those terms, the coupling between the rel¬ 
ative motion and the mass asymmetry might be related to the neck dynamics. Third, we 
used a locally harmonic approximation, but the diabatic potential might need anharmonic 
corrections. Finally, the Fokker-Planck equation can be extended to include the second 
derivatives with respect to collective coordinates (and with respect to a collective coordi¬ 
nate and the corresponding conjugate momentum) [28] . In principle, all the above effects can 
be calculated in a similar way used in this work. We hope to discuss the general analysis in 
future presentations. 
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